1 



o\ 

< 
o 

(N 



> 

O 
oo 
O 



X 



Complex zeros of the 2d Ising model on dynamical random lattices * 
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We study the zeros in the complex plane of the partition function for the Ising model coupled to 2d quantum 
gravity for complex magnetic field and for complex temperature. We compute the zeros by using the exact solution 
coming from a two matrix model and by Monte Carlo simulations of Ising spins on dynamical triangulations. We 
present evidence that the zeros form simple one-dimensional patterns in the complex plane, and that the critical 
behaviour of the system is governed by the scaling of the distribution of singularities near the critical point. 



1. THE MODEL 

We study the Ising model on dynamical square 
and triangular lattices with the spins located on 
the N faces and on the A'^ vertices, respectively. 
The partition function for a fixed lattice Gat is 
given by 
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and the partition function Z{N, (3, H) for the 
model coupled to quantum gravity is obtained 
by summing Zfiat(GAr, /9, If) over all lattices Gn 
with spherical topology. The system undergoes a 
third order phase transition 1^. Assuming scal- 
ing, the critical exponents (that are simply re- 
lated to the flat space exponents [Q ) can be com- 
puted directly from the equivalence of Z{N, (3, H) 
to the free energy of a hermitean two-matrix 
model In the usual notation, /3 = l/2,7 = 2, 
(5 = 5, vdn = 3 {dn w 4 is the Hausdorff dimen- 
sion of the system). 

If we fix (i {H), Zflat(GAr, /?, H) is a polynomial 
in y = e~^^ (c = e~^'^). In the thermodynamic 
limit its zeros (called Lee- Yang ^ or Fisher 
zeros respectively) form dense sets on lines, which 
are Stokes lines for Z. The Lee- Yang zeros lie on 
the unit circle [|l| for fixed Gn- One then defines a 
density of zeros pYL{f3,0) [H = iO) {pf{P,H) for 
Fisher zeros). In the high temperature phase the 
Lee- Yang zeros form a gap such that pYL{f3, 0) = 
for |0| < 00- P ^ f3c, the gap closes and the 
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zeros pinch the real axis at y — 1, signalling the 
onset of the phase transition. 

2. EXACT CALCULATION 

The exact solution of the Ising model on a 
dynamical lattice was given in j5|. The grand 
canonical partition function Z{g, (3, H) is given 
by Z(.g, /3, H) = Y.n=i k^Z{N, /3, H), where k = 
~4gc/{l — c^Y . Using the exact solution we de- 
termined Z{N,(3,H) for N < 14. Then its roots 
for either fixed (3 or fixed H were determined. 

Our results indicate that the Lee- Yang zeros 
lie on the unit circle (cf. [^), also after summing 
over all Gn- Fig. 1 indicates a vanishing gap 
in pyl{(3c,S) sls N oo. We also observed the 
vanishing of the gap for a fixed lattice size as /3 ^ 
Pc from the hot phase. 
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Figure 1. Lee- Yang zeros for dynamical lattices 
of varying size N = 8, 9, 14 at /3c- 

Also the Fisher zeros form Id patterns. For 

H = and N ^ oo they approach the (physical 
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and non-physical) critical points c = ±1/4. There 
are also Fisher zeros on the imaginary axis, that 
flow to c = ±ioo. Fig. 2 shows the Fisher zeros 
for H = and H at fixed A''. Mapping the ze- 
ros onto the c-plane, where the tilde refers to the 
dual spin model (recall the usual duality relation 
c = e~'^^ = tanh(/?)), we observed the existence 
of an antiferromagnetic phase transition in the 
dual model in agreement with ref. This AF 
transition is absent in the original model where 
the spins are located on the faces of the lattice. 



(a) c-plane 



(b) tanh(p)-plane 
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Figure 2(a). Fisher zeros in the complex c plane 
for fixed = 14 and varying y. Only for y = 
1 {H = 0) the zeros pinch the real axis, (b) 
The trajectories in the c-plane. Note the FM and 
AFM phase transitions at c = 3/5, 5/3 for y — I. 
The approach to these points is observed as N 
increases (not shown). 

A priori, it is not evident that the partition 
function zeros will exhibit scaling in the presence 
of quantum gravity. However, the evidence in the 
range of central charge < c < 1 for a diverg- 
ing correlation length at the critical point makes 
such a scenario more plausible. 

From scaling arguments ||] one deduces that 
the jth Lee- Yang zero should scale as Hj ^ 
j^-f35/(vdH)_ Yig. 3 shows the first three Lee- 
Yang zeros vs. N. Fig. 3 gives l35/{vdH) ~ 
0.871 ± 0.002, 0.935 ± 0.002 and 0.951 ± 0.002 
for j = 1,2,3 respectively. The value expected 
from the scaling assumption is 5/6 w 0.8333... for 
ah J. 

It is expected from renormalization group ar- 
guments ||] that the trajectories of the motion 
of Fisher zeros in the K plane [K = u ~ Uc and 
u = e~'^^) with varying real H, will form an an- 
gle i(; — tt/ (2/3(5) with the real K axis as ^ oo. 
For our case ip — 36°. The approach to this value 
(dashed line) can be observed in Fig. 4. 



Next we tested the scaling relation g 
i?|(Af/j)2*/(*+i' = F (if(7V/jy/('^'^")) where F 
is a universal scaling function. Taking K real, 
we plotted Hj{N/jf/^ vs. {(3 - l3c){N / jf/^ (not 
shown) for all zeros of order j > 1 and all lattice 
sizes. The points fell on one universal curve rep- 
resenting the function F for large j in the hot 
phase. By varying the value of the exponent 
5/6 we observed a broadening of the curve and 
obtained a (somewhat subjective) determination 
S/{5 -M) = (36/{vdH) ^ 0.85 ± 0.05, in excellent 
agreement with the KPZ value. 

From the definition of F we expect for the 
jth Fisher zero and large iV: InjiiTjl ~ ln|/3j — 
(3c\ = l/ivdH)Hj/N) + C Fig. 5 shows \Ki\ 
and \Pi — /3c I vs. N on a log-log plot. Kj did 
not yield a straight line at these lattice sizes, but 
\Pj-f3c\ scales and gives l/(i^di/) « 0.327±0.001, 
0.377 ± 0.002, for j = 1, 2 respectively The value 
expected from scaling is 1/3. 
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Figure 3. Scaling of Lee- Yang zeros Hj 
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Figure 4. tantp versus N. The dashed line is 
tan36°. 
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Figure 5. Scaling of the first Fisher zero. 
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3. MULTIHISTOGRAMMING 

For systems with more than a few spins, the 
study of complex zeros must rely on numeri- 
cal methods. We used multihistogramming [|j 
to determine the partition function accurately 
throughout a continuous range of couplings. The 
lattice sizes that we simulated ranged from 32 
to 256 vertices (60 - 508 triangles). We refer to 
|lO[ for details. Lee- Yang and Fisher zeros were 
determined from the minima of Z {(3 , H) / Z {(3) 
and Z{(3,0)l Z{Ke(3), respectively, for real [3 and 
imaginary H and for complex fj and H — 0. 
Tables with the observed zeros can be found in 
|l0| . We observed Lee- Yang zeros as follows: for 
Ny — 64, only the first zero could be seen, for 
Nv = 96, the first two zeros, for = 128, three 
zeros, for Ny = 256, five zeros. For the Fisher 
zeros, only the first zero was visible for all lattice 
sizes Ny = 32, 64, 96, 128, 256. 
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Figure 6. The first and second Lee- Yang zero 
observed using multihistogramming. 
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Figure 7. The first Fisher zero observed using 
multihistogramming. 

Fits to Hj ~ N-PS/{^d„) gj^g foj. ^Yie first 
and second Lee- Yang zero (Fig. 6) l35/{vdH) ~ 
0.773 ±0.013, 0.788 ±0.033 (j = 1,2). Using 7V„ 
instead of N for volume, the corresponding values 
are0.787±0.013, 0.800±0.034 (j = 1,2), indicat- 
ing finite size effects. The results are in quite good 
agreement with the expected value 5/6. From 



z^dff = 3 § we obtain /3(5 = 2.36±0.04, 2.40±0.10 
[j — 1,2 respectively); the exact value is 2.5. 

The scaling of the observed Fisher zeros is 
shown in Fig. 7. The extracted values of l/{vdH) 
are: for K 0.392 ± 0.016; for |/3 - /3c| 0.386 ± 

0. 014 (expected value: 1/3 for both as TV ^ cxo). 
They increase slightly if Ny is used instead of N . 
Including only the three largest lattices we obtain 
(using N for volume): l/iydn) = 0.350 ± 0.067. 

4. CONCLUSION 

We conclude that the critical behaviour of the 
system is described by the scaling of the dis- 
tribution of the complex zeros of the partition 
function. Our result for Lee- Yang zeros presents 
us with the challenge of proving a corresponding 
Lee- Yang theorem for the case when a fluctuating 
metric contributes an additional quantum degree 
of freedom. 
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